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1 Heavy Electrons 

1.1 Introduction 

In a world where it is possible to hold a levitated high temperature supereonduetor in the palm of 
one’s hand, it is easy to forget the ongoing importanee of low temperature researeh. Heavy elec¬ 
tron materials are a class of strongly correlated electron material containing localized magnetic 
moments which, by entangling with the surrounding electrons, profoundly transform the metal¬ 
lic properties. A heavy fermion metal can develop electron masses 1000 times that of copper, it 
can also develop unconventional superconductivity, transform into new forms of quantum order, 
exhibit quantum critical and topological behavior. Although most of these properties develop 
well below the boiling point of nitrogen, the diversity and highly tunable nature of their ground- 
states make them an invaluable vital work-horse for exploring and researching the emergent 
properties of correlated quantum matter. 

This lecture will give an introduction to heavy fermion materials, trying to emphasize a 21st 
century perspective. More extensive discussion and development of the ideas in these notes 
can be found in an earlier review article [|T|| and the latter chapters of my book “Introduction to 
Many Body Physics” Q. 


Increasing localization 


Magnetic moments 


4f 

Ce 

Pr 

Nd 

Pm 

Sm Eu 

6d 

Tb^ 




5f 

Th 

Pa 

U 

Np 

Pu Am 

Cm Bk Cf 

Es 

Fm ^ 


3d 

Sc 


Ti 

V 

Cr 

Mn 

Fe 

Co 

Ni 

Cu 

4d 

Y 


Zr 

Nb 

Mo 

Tc 

Ru 

Rh 

Pd 

Ag 

5d 

Lu 


HI 

Ta 

W 

Re 

Os 

Ir 

Pt 

Au 













increasing 

localization 


Superconductivity 



transition actinide rare earth 


Fig. 1: The Kmetko-Smith diagram showing the broad trends towards increasing electron 
localization in the d- and f-electron compounds. 

In the periodic table, the most strongly interacting electrons reside in orbitals that are well- 
localized. In order of increasing localization, partially filled orbitals are ordered as follows: 

5d <4d <3d <5f < 4/. (1) 

In addition, when moving along a row of the periodic table, the increasing nuclear charge 
pulls the orbitals towards the nucleus. These trends are summarized in the “Kmetko-Smith 
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diagram” Q in Fig The d-orbital metals at the bottom left of this diagram are highly itiner¬ 
ant and exhibit conventional superconductivity. By contrast, in rare earth and actinide metals 
towards the top right-hand comer, the f-shell electrons are localized, forming magnets or anti- 
ferromagnets. It is the materials that lie in the cross-over between these two regions that are 
particularly interesting, for these materials are “on the brink of magnetism”. It is in this cross¬ 
over region that many strongly correlated materials reside: it is here for instance, that we find 
cerium and uranium, which are key atoms for a wide range of 4f and 5f heavy electron materials. 

1.2 Local moments and the Kondo effect 

Heavy electron materials contain a lattice of localized electrons immersed in a sea of mobile 
conduction electrons. To understand their physics, we need to first step back and discuss in¬ 
dividual localized moments, and the mechanism by which they interact with the surrounding 
conduction sea. 

The key feature of a localized moment, is that the Coulomb interaction has eliminated the high 
frequency charge fluctuations, leaving behind a low energy manifold of degenerate spin states. 
In rare earth and actinide ions, the orbital and spin angular momentum combine into a single 
entity with angular momentum j = I + s. For example, a ion contains a single unpaired 
4f-electron in the state 4/\ with I = 3 and 5 = 1/2. Spin-orbit coupling gives rise to low- 
lying multiplet with j = 3 - \ = |, consisting oflj + \ =6 degenerate orbitals |4/^ : Jm), 
(nij 6 [-|, |]) with an associated magnetic moment M = 2.64/^^. In a crystal, the 2j + I fold 
degeneracy of such a magnetic ion is split, and provided there are an odd number of electrons 
in the ion, Kramer’s theorem guarantees that the lowest lying state has at least, a two fold 
degeneracy. (Fig. |^a and b.) 

One of the classic signatures of localized moments, is a high temperature Curie Weiss suscepti¬ 
bility, given by 

^ * ^‘ 3(7 + 6 ) + 1 ). ( 2 ) 

where, n,- is the concentration of magnetic moments while M is the magnetic moment with total 
angular momentum quantum number j and gyro-magnetic ratio (“g-factor”) g. 6 is the “Curie 
Weiss” temperature, a phenomenological scale which takes account of interactions between 
spins. 

The presence of such local moments inside a metal profoundly alters its properties. The physics 
of an isolated magnetic ion is described by the Kondo model 

AH 

^ ........ 

H = • ( 3 ) 

ko- 

where c\a- creates a conduction electron of energy ek, momentum k and i//\0) = <^\cr 

creates a conduction at the origin, where Ns is the number of sites in the lattice. The conduction 
sea interacts with local moment via an antiferromagnetic contact interaction of strength J. The 
antiferromagnetic sign (7 > 0) of this interaction is an example of “super-exchange”, first 
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Fig. 2: (a) In isolation, the localized atomic states of an atom form a stable, sharp excitation 
lying below the continuum, (b) In a crystal, the 2j+\ fold degenerate state splits into multiplets, 
typically forming a low lying Kramers doublet, (c) The inverse of the Curie-Weiss susceptibility 
of local moments x~^ Is a linear function of temperature, intersecting zero at T = -9. 


predieted by Philip W. Anderson whieh results from high energy valenee fluetuations. 

Jun Kondo Q first analyzed the effeet of this seattering, showing that as the temperature is 
lowered, the effeetive strength of the interaetion grows logarithmieally, aeeording to 

J ^ J{T) = J + 2J^p\n^ (4) 

where p is the density of states of the eonduetion sea (per spin) and D is the band-width. The 
growth of this interaetion enabled Kondo to understand why in many metals at low temperatures, 
the resistanee starts to rise as the temperature is lowered, giving rise to resistance minimum. 



Fig. 3: (a) Schematic temperature-field phase diagram of the Kondo effect. At fields and tem¬ 
peratures large compared with the Kondo temperature Tk, the local moment is unscreened with 
a Curie susceptibility. At temperatures and fields small compared with Tk, the local moment is 
screened, forming an elastic scattering center within a Landau Fermi liquid with a Pauli sus¬ 
ceptibility x ~ 7 ^- (b) Schematic susceptibility curve for the Kondo effect, showing cross-over 
from Curie susceptibility at high temperatures to Pauli susceptibility at temperatures below the 
Kondo temperature Tk- (c) Specific heat curve for the Kondo effect. Since the total area is the 
full spin entropy R\n2 and the width is of order T^, the height must be of order y ~ R\n2ITx- 
This sets the scale for the zero temperature specific heat coefficient. 

Today, we understand this logarithmie eorreetion as a renormalization of the Kondo eoupling 
constant, resulting from fact that as the temperature is lowered, more and more high frequency 
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quantum spin fluetuations beeome eoherent, and these strengthen the Kondo interaetion. The 
effeet is closely analogous to the growth of the strong-interaction between quarks, and like 
quarks, the local moment in the Kondo effect is asymptotically free at high energies. However, 
as you can see from the above equation, once the temperature becomes of order 


Tk 


Dexp 


1 


the correction becomes as large as the original perturbation, and at lower temperatures, the 
Kondo interaction can no longer be treated perturbatively. In fact, non-perturbative methods 
tell us that this interaction scales to strong coupling at low energies, causing electrons in the 
conduction sea to magnetically screen the local moment to form an inert Kondo singlet denoted 
by 

|GS) = ^(IM)-|llT)), (5) 

V2 

where the thick arrow refers to the spin state of the local moment and the thin arrow refers to the 
spin state of a bound-electron at the site of the local moment. The key features of the impurity 
Kondo effect are are: 


• The electron fluid surrounding the Kondo singlet forms a Fermi liquid, with a Pauli sus¬ 
ceptibility x ^/Tk- 


• The local moment is a kind of qubit which entangles with the conduction sea to form 
a singlet. As the temperature T is raised, the entanglement entropy converts to thermal 
entropy, given by the integral of the specific heat coefficient. 


SiT) = 



dT 


, Cv{T') 

r 


Since the total area under the curve, S{T ^ oo) = i? In 2 per mole is the high tempera¬ 
ture spin entropy, and since the characteristic width is the Kondo temperature, it follows 
that the the characteristic zero temperature specific heat coefficient must be of order the 


inverse Kondo temperature: y = ^(T 0) 


R\n2 

Tk 


. (See Fig. 3 b) 


• The only scale in the physics is Tk- For example, the resistivity created by magnetic 
scattering off the impurity has a universal temperature dependence 


RjT) 

Ru 



( 6 ) 


where n, is the concentration of magnetic impurities, <P(x) is a universal function and pu 
is the unit of unitary resistance (basically resistance with a scattering rate of order the 
Fermi energy). 


Ri 


2ne^ 

nmp 


(7) 


Experiment confirms that the resistivity in the Kondo effect can indeed be scaled onto a 
single curve that fits forms derived from the Kondo model (see Fig|^. 











5.6 


Piers Coleman 


• The seattering off the Kondo singlet is resonantly eonfined to a narrow region of order 
Tk, ealled the Kondo or Abriksov-Suhl resonance. 
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Fig. 4: Temperature dependence of resistivity associated with scattering from an impurity 
spin from m- The resistivity saturates at the unitarity limit at low temperatures, due to the 
formation of the Kondo resonance. Adapted from 


1.3 The Kondo lattice 

In heavy fermion material, containing a lattice of local moments, the Kondo effect develops 
coherence. In a single impurity, a Kondo singlet scatters electrons without conserving momen¬ 
tum, giving rise to a huge build-up of resistivity at low temperatures. However, in a lattice, with 
translational symmetry, this same elastic scattering now conserves momentum, and this leads to 
coherent scattering off the Kondo singlets. In the simplest heavy fermion metals, this leads to a 
dramatic reduction in the resistivity at temperatures below the Kondo temperature. 

As a simple example, consider CeCug a classic heavy fermion metal. Naively, CeCug is just 
a copper alloy, in which 14% of the copper atoms are replaced by cerium, yet this modest 
replacement radically alters the metal. In this material, it actually proves possible to follow the 
development of coherence from the dilute single ion Kondo limit, to the dense Kondo lattice, by 
forming the alloy Tdi\^fZeyCu^. Lanthanum is iso-electronic to cerium, but has an empty f-shell, 
so the limit jc ^ 0 corresponds to the dilute Kondo limit, and in this limit the resistivity follows 
the classic Kondo curve. However, as the concentration of cerium increases, the resistivity 
curve starts to develop a coherence maximum, an in the concentrated limit drops to zero with a 
characteristic T^ dependence of a Landau Fermi liquid (see Fig. |^. 

CeCug displays the following classic features of a heavy fermion metal: 

• A Curie-Weiss susceptibility ~ {T + Q)~^ at high temperatures. 

• A paramagnetic spin susceptibility^ ~ cons at low temperatures. 
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• A dramatieally enhaneed linear speeifie heat Cy = jT at low temperatures, where in 
CeCue y ~ lOOOmJ/mol/K^ is about 1000 times larger than in eopper. 

• A quadratie temperature dependenee of the low temperature resistivity p = p,, + AT^ 


In a Landau Fermi liquid Q, the magnetie susceptibility and the linear specific heat coefficient 
7 = Cy/r|r_o are given by 


A 

7 


(Msf 


N*(0) 


3 


N*{0) 


( 8 ) 


where N*{0) = ^N(0) is the renormalized density of states and is the spin-dependent part 
of the s-wave interaction between quasiparticles. One of the consequences of Fermi liquid 
theory, is that the density of states factors out of the Sommerfeld or Wilson ratio between the 
susceptibilty and linear specific heat coefficient. 


7 \2nkBl l+F’" 


(9) 


In heavy fermion metals, this ratio remains approximately fixed across several decades of vari¬ 
ation in X and y. This allows us to understand heavy fermion metals as a lattice version of the 
Kondo effect gives rise to a renormalized density of states N*{0) ~ 

The discovery of heavy electron compounds in the 1970s led Mott |10| and Doniach pT| to 
propose that heavy electron systems should be modeled as a “Kondo-lattice”, where a dense 
array of local moments interact with the conduction sea via an antiferromagnetic interaction 
J. In such a lattice, the local moments polarize the conduction sea, and the resulting Friedel 
oscillations in the magnetization give rise to an antiferromagnetic RKKY (Rudermann Kittel 
Kasuya Yosida) magnetic interaction p2] - [T4| that tends to order the local moments. Mott and 
Doniach realized that this interaction must compete with the Kondo effect. 

The simplest Kondo lattice Hamiltonian [ 15 |is 

H — ^ CkC kcr^-kiT + J ^ ^ j ‘ ^ jo^afiCjp, (10) 


k(T 


where 


■ - 
^ ja — 




( 11 ) 


creates an electron at site j. Mott and Doniach p0|[TT] | pointed out that there are two en¬ 
ergy scales in the Kondo lattice: the Kondo temperature Tk ~ and the RKKY scale 

Erkky = Erkky = J^P- For small Jp, Erkky » Tr leading to an antiferromagnetic ground- 
state, but when Jp is large, Tr » Errry, stabilizing a ground-state in which every site in the 
lattice resonantly scatters electrons. Based on a a simplified one-dimensional “Kondo neck¬ 
lace” model [ 16|, Doniach conjectured pT| that the transition between the antiferromagnet and 
the dense Kondo ground state is a continuous quantum phase transition. Experiment confirms 
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Fig. 5: Doniach phase diagram for the Kondo lattice, illustrating the antiferromagnetic regime 
and the heavy fermion regime, for Tk < Trkky cmd Tk > Trkky respectively. The effective 
Fermi temperature of the heavy Fermi liquid is indicated as a solid line. Experimental evidence 
suggests that in many heavy fermion materials this scale drops to zero at the antiferromagnetic 
quantum critical point. 


this eonjeeture, and today we have several examples of sueh quantum eritieal points, ineluding 
CeCug doped with gold to form CeCug-j^Aujc and CeRhIns under pressure [ 17 -191. In the fully 
developed Kondo lattiee ground state Bloeh’s theorem insures that the resonant elastie seat- 
tering at eaeh site will generate a renormalized f- band, of width ~ Tk- In eontrast with the 
impurity Kondo effect, here elastic scattering at each site acts coherently. For this reason, as the 
heavy electron metal develops at low temperatures, its resistivity drops towards zero (see Fig. 

ig». 


In a Kondo lattice, spin entanglement is occurring on a truly macroscopic scale, but this entan¬ 
glement need not necessarily lead to a Fermi liquid. Experimentally, many other possibilities 
are possible. Here are some examples. 


• CcsBuPts, a Kondo insulator in which the formation of Kondo singlets with the Ce mo¬ 
ments drives the development of a small insulating gap at low temperatures and 

• CeRhInj, an antiferromagnet on the brink of forming a Kondo lattice, which under pres¬ 
sure becomes a heavy fermion superconductor with rc=2K. 

• UBcia a heavy fermion superconductor which transitions directly from an incoherent 
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metal with resistivity 200/ii3em, into a supereondueting state. 

Eaeh of these materials has qualitatively the same same high temperature Curie Weiss mag¬ 
netism and the same Kondo resistivity at high temperatures, due to ineoherent scattering off the 
local moments. However at low temperatures the scattering off the magnetic Ce ions becomes 
coherent and new properties develop. 

^ 100 200 




Fig. 6: (a) Resistivity of CexLa\-xCu^. Dilute Ce atoms in LaCu^ exhibit a classic “Kondo” 
resistivity, but as the Ce concentration becomes dense, elastic scattering off each Ce atom leads 
to the development of a coherent heavy fermion metal, (b) Resistivities of four heavy fermion 
materials showing the development of coherence. A variety of antiferromagnetic magnetic, 
Fermi liquid, superconducting and insulating states are formed (see text). 


2 Kondo insulators: the simplest heavy fermions 


In many ways, the Kondo insulator is the simplest ground-state of the Kondo lattice. The 


first Kondo insulator (KI), SmBgwas discovered almost fifty years ago [20| and today there 
are several known examples including CeaBUPta. At room temperature, these KIs are metals 
containing a dense array of magnetic moments, yet on cooling they develop a narrow gap due 
the formation of Kondo singlets which screen the local moments [pT]-(2^. We can gain a lot of 


insight by examining the strong coupling limit in which the dispersion of the conduction sea is 
much smaller than the Kondo coupling J. Consider a simple tight-binding Kondo lattice 

= (c jpCr ffixCja) 


CTi 


( 12 ) 


H = -t^ {c\a-Cjo- + H.C) + § j, 

j,a/3 

in which t/J << 1 is a small parameter. In this limit, the inter-site hopping is a perturbation to 
the on-site Kondo interaction, 

//y-^o 


H 




(13) 
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and the eorresponding ground-state corresponds to the formation of a spin singlet at each site, 
denoted by the wavefunction 

\Ki) = [] (14) 

where the double and single arrows denote the localized moment and conduction electron re¬ 
spectively. 

Each singlet has a ground-state energy E = -\j per site and a singlet-triplet spin gap of magni¬ 
tude AE = 27. Moreover, if we either remove an electron from site i, we break a Kondo singlet 
and create an unpaired spin with excited energy |7, 

lqp^ i T) n = V2c,i|7:/), (15) 

as illustrated in Fig|7];a). Similarly, if we add an electron, we create an electron quasiparticle, 
corresponding to an unpaired local moment and a doubly occupied conduction electron orbital 

lqp“, i T) (r/ii) n - U;T;) = (16) 

j*i 

as illustrated in Fig|7];b). 


Hole doping 

Electron doping 



Fig. 7: Showing (a) hole and (b) electron doping of strong coupling Kondo insulator, (c) 
Dispersion of strong coupling Kondo insulator. A small amount of hold doping 6 gives rise to a 
“large ” Eermi surface containing 2 — 6 heavy electrons. 

If we now reintroduce the hopping -t between sites, then these quasiparticle excitations become 
mobile, as illustrated in Fig. [TJ^a) and (b). From the explicit form of the states, we find that the 
nearest neighbor hopping matrix elements are (qp*, zcr|//|qp*, ycr) = +|„ giving quasiparticle 
energies 


3 


(17) 
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To transform from the quasipartiele, to the eleetron basis, we need to reverse the sign of the 
hole (qp"^) dispersion to obtain the valence band dispersion, so that the band energies predicted 
by the strong coupling limit of the Kondo lattice are 


~ ti^Cx + Cy + C^) + ^J, 


(18) 


separated by an energy 37 as shown in Fig. |^c). Note that these are “hard core” fermions that 
can not occupy the same lattice site simultaneously. 

In this way, the half-filled strong coupling Kondo lattice forms an insulator with a charge gap 
of size 37 and a spin gap of size 27. Notice finally that if we dope the insulator with an amount 
6 of holes, we form a band of heavy fermions. In this way, Kondo insulators can be considered 
the parent states of heavy electron materials. However, we’d like to examine the physics of a 
Kondo lattice at weak coupling, and to do this requires a different approach. 


3 Large N Expansion for the Kondo Lattice 


3.1 Philosophy and Formulation 

One of the great difficulties with the Kondo lattice, is that there is no natural small parame¬ 
ter to carry out an approximate treatment. One way around this difficulty, is to use a large N 
expansion, in which we extend the number of spin components of the electrons from 2 to N. 
Historically, Anderson [j^ pointed out that the large spin-orbit coupling in heavy fermion com¬ 
pounds generates (if we ignore crystal fields) a large spin degeneracy N = 2j + I, furnishing a 
small parameter 1 /N for a controlled expansion about the limit N ^ oo. One of the observations 


arising from Anderson’s idea [26 27 j is that the RKKY interaction becomes negligible (of or¬ 
der 0(1/N^)) in this limit and the Kondo lattice ground-state becomes stable. This observation 


opened the way to path integral mean-field treatments of the Kondo lattice [27-32| 

The basic idea of the large N limit is to examine a limit where every term in the Hamiltonian 
grows extensively with N. In the path integral for the partition function, the corresponding 
action then grows extensively with N, so that 


= J !D[ilj]e = J' 2)i/^exp = J £)[i/r] exp 




eff _ 


(19) 


Here ^ ~ figff behaves as an effective Planck’s constant for the theory, focusing the path integral 
into a non-trivial “semi-classical” or “mean field” solution as fieff ^0. As A ^ oo, the 
quantum fluctuations of intensive variables a, such as the electron density per spin, become 
smaller and smaller, scaling as {6a^)/{a^) ~ 1/A, causing the path integral to focus around a 
non-trivial mean-field trajectory. In this way, one can obtain new results by expanding around 
the solvable large A limit in powers of (Fig. 


We will use a simplified Kondo lattice model introduced by Read and Newns [27 j, in which all 
electrons have a spin degeneracy A = 2j -l- 1, 

H ^ ko-^kcy ^ 'j C ja^ afid) • 


kof 


7,0-/3 


( 20 ) 
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Fig. 8: Illustration of the convergence of a quantum path integral about a semi-classical trajec¬ 
tory in the large N limit. 

where ereates an eleetron loealized at site j and the spin of the loeal 

moment at position Ry is represented by pseudo-fermions 

Sa/^U)=fjah-^f-Sa^. ( 21 ) 

This representation requires that we set a value for the eonserved / oeeupaney nf(j) = Q at 
each site. This interaction can be rewritten in a factorized form 


ktT 


t _l. 

ko'^kci' 


^ {f^ • 

.hap 


( 22 ) 


Read Newns model for the Kondo lattice 


where the potential scattering terms resulting from the rearrangement of the f-operators have 
been absorbed into a shift of the chemical potential. Note that: 


the model has a global S UfN) symmetry associated with the conserved magnetization. 


the Read Newns (RN) model is a lattice version of the Coqblin-Schrieffer Hamiltonian 


1331 introduced to describe the Kondo interaction in strongly spin-orbit coupled rare-earth 


ions. While the Coqblin-Schrieffer interaction is correct at each site, the assumption that 
the S U(N) spin is conserved by electron hopping is an oversimplification. (This is a price 
one pays for a solvable model.) 


in this factorized form, the antiferromagnetic Kondo interaction is “attractive”. 
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• the coupling constant has been scaled to vary as JjN, to ensure that the interaction grows 
extensively with N. The interaction involves a product of two terms that scale as 0{N), 
so that J/N X 0{N^) ~ 0{N). 

• the RN model also has a local gauge invariance: the absence of f-charge fluctuations 
allows us to change the phase of the f-electrons independently at each site 

fj. ^ (23) 


A tricky issue concerns the value we give to the conserved charge nf = Q. In the physical 
models of interest, nf = 1 at each site, so one might be inclined to explicitly maintain this 
condition. However, the large N expansion requires that the action is extensive in N, and this 
forces us to consider more general classes of solution where Q scales with N so that the filling 
factor q = QjN is finite as N ^ oo. Thus if we’re interested in a Kramer’s doublet Kondo 
model, we take the half-filled case q = lj2, Q = N/2, but if we want to understand a 7 = 7/2 
Yb^'^ atom without crystal fields, then in the physical system A = 27 -l -1 = 8 , and we should fix 
q = Q/N=l/S. 

The partition function for the Kondo lattice is then 


Z = 


Tr 


e PI 6(hf( j) - Q) 

j 


(24) 


where 6(hf(j) - Q) projects out the states with nf{i) = Q at site 7 . By re-writing the delta 
function as a Fourier transform, the partition function can be can be rewritten as a path-integral. 



L[i//',i//,A] 

j" d]exp 

i H ^ 

- 1 drif/-t-H\ijj,il/,Ai^ 

^0 


(25) 


where = ({c^}, {/^}) schematically represent the conduction and f-electron fields, 

= x; kor^kti' 

ka j,a/3 j 

The field Aj is a fluctuating Lagrange multiplier that enforces the constraint nj = Q at each site. 
Next we carry out a Hubbard Stratonovich transformation on the interaction, 

“ ^ JAfjA) ^ ^ \^j jafja) + jaCja) Tyj -|- N—j—. (27) 

ap a 

In the original Kondo model, we started out with an interaction between electrons and spins. 
Now, by carrying out the Hubbard Stratonovich transformation, we have formulated the inter¬ 
action as the exchange of a charged boson 
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-i Z (28) 

where the solid lines represent the conduction electron propagators, and the dashed lines repre¬ 
sent the f-electron operators. Notice how the bare amplitude associated with the exchange boson 
is frequency independent, i.e the interaction is instantaneous. Physically, we may interpret this 
exchange process as due an intermediate valence fluctuation. 

The path integral now involves an additional integration over the hybridization fields V and V, 


S[V,V,A 


J D[y,y,T] J D[iA^<A]exp[-j^V5r<A + ^[v,y,d]) 


H[V, y T] = 6kC^k<rCk. + Xi Vj + Aj(nfj - Q) + Ay^|?9) 


J 


Read Newns path integral for the Kondo lattice. 


where we have suppressed summation signs for repeated spin indices (summation convention). 
The RN path integral allows us to develop a mean-field description of the many body Kondo 
scattering processes that captures the physics and is asymptotically exact as N ^ oo. In this 
approach, the condensation of the hybridization field describes the formation of bound-states 
between spins and electrons that can not be dealt with in perturbation theory. Bound-states 
induce long range temporal correlations in scattering: once the hybridization condenses, the 
interaction lines break-up into independent anomalous scattering events, denoted by 



The hybridization V in the RN action carries the local [/(I) gauge charge of the f-electrons, 
giving rise to an important local gauge invariance: 


fj<r ^ Vj ^ Aj ^ Aj - (30) 

Read Newns gauge transformation. 


This invariance can be used to choose a gauge in which Vj is real, by absorbing the phase of 
the hybridization Vj = \Vj\e“^j into the f-electron. In the radial gauge. 
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s[|y|4 


Z = J mv\,A] J D[iA^<A]exp[- + 


|V;f 


(31) 


Read Newns path integral: “radial gauge”. 


Subsequently, when we use the radial gauge, we will drop the moduli! sign. The interesting 
feature about this Hamiltonian, is that with the real hybridization, the eonduetion and f-eleetrons 
now transform under a single global f/(l) gauge transformation, i.e the f-electrons have become 
charged. 


3.2 Mean-Field Theory 


The interior fermion integral in the path integral ( [3T] ) defines an effective action SeW, d] by the 
relation 


ZE = exp[-NSE[V,A]] = J ©[<A^<A]exp[-5[y,d, , 


(32) 


The extensive growth of the effective action with N means that at large N, the integration in (29) 
is dominated by its stationary points, allowing us to dispense with the integrals over V and A. 


-f 


n[A, V] exp [-NSeIV, d]] « exp [-iV5£[y, d]] 


Saddle Point 


(33) 


In practice, we seek uniform, static solutions, y/r) = y Aj(T) = A. In this case the saddle 
point partition function Ze = Tre^d^^wFr simply the partition function of the static mean-field 
Hamiltonian 


h(k) 


xfek 

y' 

^kcr 

l[v 

d. 



HmFT - ^ f^kcr) ^ 

ktr ^ 

— ^ ^ k(T ^(k) l/^kc 


kcr 


+ m|^-d<5rj 
+ m|^-d<5rj. 


(34) 


Here, f\a- = Z;is the Fourier transform of the /-electron field and we have 
introduced the two component notation 


/ka 


^kcr 


.fk< 


\Jkcr/ 


/ kcr — ko-, / ,k(T ^ ? ^(k) 


^6k 

y d 


(35) 


We should think about Hmft as a renormali z ed Hamiltonian, describing the low energy quasi¬ 
particles, moving through a self-consistently determined array of resonant scattering centers. 
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Fig. 9: (a) Dispersion for the Kondo lattice mean field theory, (b) Renormalized density of 
states, showing “hybridization gap” (Ag). 


Later, we will see that the f-eleetron operators are composite objects, formed as bound-states 
between spins and conduction electrons. 

The mean-field Hamiltonian can be diagonalized in the form 


HmFT — ^ ^k(T; 

kcr 


^kcr 


0 


J 


+ NK 


,"kcr J 


vv 

~T 


Aq . 


(36) 


Here a\cr = + Vk/\o- and b\cr = -VkC^ko- + Hkf\a- are linear combinations of c^ko- and 

f\a; playing the role of “quasiparticle operators” with corresponding energy eigenvalues 


Det 

^kl - 

S_ 



(y Aj\ 


= (Fk± - 6k)(Fk± - A)-\V\^ = 0, 


or 


Ek± 


fk + d 


fk - d 


+ |V|^ 


and eigenvectors taking the BCS form 


j Mk [ _ 

1 , (6k-d)/2 

1 J 

^ 2 ^(2^)' + |y|2 


The hybridized dispersion described by these energies is shown in Fig. 
Note that: 


(37) 


(38) 


(39) 


• The Kondo effect injects an f-band into the conduction sea, hybridizing with the conduc¬ 
tion band to create two bands separated by a direct “hybridization gap” of size 2V and a 
much smaller indirect gap. If we put Ck = ±7), we see that the upper and lower edges of 
the gap are given by 


+D +A 


'+D-A 


+ y 2 




(D » A) 


(40) 
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so the indireet gap has a size Ag ~ 2V^/D, where D is the half-bandwidth. We will 
see shortly that V^/D ~ T/f is basieally the single-ion Kondo temperature, so that that 
V ~ is the geometrie mean of the band-width and Kondo temperature. 

• In the case when the chemical potential lies in the gap, a Kondo insulator is formed. 

• A conduction sea of electrons has been transformed into a heavy Fermi sea of holes. 

• The Fermi surface volume expands in response to the formation of heavy electrons (see 
Fig. [T^ to count the total number of occupied quasiparticle states 

Ntot = nuAa) = («/ + flc) (41) 

kA(T 

where nua- = (^kAa-akAa- is the number operator for the quasiparticles and ric is the total 
number of conduction electrons. This means 

Ntot = N-^ = Q + n,, (42) 

(Iny 

where is the volume of the unit cell. This is rather remarkable, for the expansion of the 
Fermi surface implies an increased negative charge density in the Fermi sea. Since charge 
is conserved, we are forced to conclude there is a compensating +Q\e\ charge density per 
unit cell provided by the Kondo singlets formed at each site, as illustrated in Fig. [T^ 



Fig. 10: (a) High temperature state: small Fermi surface with a background of spins; (b)Low 
temperature state where large Fermi surface develops against a background of positive charge. 
Each spin “ionizes” into Q heavy electrons, leaving behind a background of Kondo singlets, 
each with charge +Qe. 


3.3 Free energy and Saddle Point 


Let us now use the results of the last section to calculate the mean-field free energy Fmft and 
determine, self-consistently the parameters A and V which set the scales of the Kondo lattice. 
By diagonalizing the mean field Hamiltonian, we obtain 


F 

N 



+ K 



(43) 
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Let us diseuss the ground-state, in whieh only the lower-band eontributes to the Free energy. As 
T ^ 0, we can replace -T ln(l ^ 0(-£’k)£’k, so the ground-state energy Eq = F{T = 0) 

involves an integral over the occupied states of the lower band: 

(44) 


where we have introduced the density of heavy electron states p*{E) = 2k,± 5{E - Now 

by ( [T7| ) the relationship between the energy E of the heavy electrons and the energy e of the 
conduction electrons is 


As we sum over momenta k within a given energy shell, there is a one-to-one correspondence 
between each conduction electron state and each quasiparticle state, so we can write p*(E)dE = 
p{e)de, where the density of heavy electron states 


p\E) 




y 2 \ 


(45) 


Here we have approximated the underlying conduction electron density of states by a constant 
p = 1/(2D). The originally flat conduction electron density of states is now replaced by a 
“hybridization gap”, flanked by two sharp peaks of width approximately npV^ ~ Tk (Fig. |^. 
Note that the lower band-width is lowered by an amount -V^jD. With this information, we can 
carry out the integral over the energies, to obtain 


(46) 

NN, I (E-AfI V 

where we have assumed that the upper band is empty, and the lower band is partially filled. 
Carrying out the integral we obtain 



where we have replaced d = npV^, which is the width of an isolated f-resonance, and have 
dropped terms of order 0{A^ jD). We can rearrange this expression, absorbing the band-width 

_ J_ 

D and Kondo coupling constant into a single Kondo temperature Tk = De ■'p as follows 
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n 



(48) 


This describes the energy of a family of Kondo lattice models with different J(D) and cutoff 
D, but fixed Kondo temperature. If we impose the constraint ^ = (uf) - Q = 0 we obtain 
^ ^ = 0, so 


EoiV) A. 

-= — In 

NN, n 


A \ D^p 


(A = np\V\^) 


(49) 


[n-qeTx) 2 

Let us pause for a moment to consider this energy functional qualitatively. There are two points 


E[V] 



Fig. 11: Mexican hat potential for the Kondo Lattice, evaluated at constant (nf) = Q as a 
function of a complex hybridization V = \ V\e’‘^ 


to be made 


The energy surface Eo(V) is actually independent of the phase of V = |y|e''^(see Fig. [IT] ), 
and has the form of “Mexican Hat” at low temperatures. The minimum of this functional 
will then determine a family of saddle point values V = \Vo\e“^, where f can have any 
value. If we differentiate the ground-state energy with respect to zl, we obtain 


0 = - In 
n 


A 


nqlK 


or 


A = nqTK 


confirming thatzl ~ Tk- 


• The mean-field value of the constraint field A is determined relative to the Fermi energy p. 
Were we to introduce a slowly varying external potential field to the conduction electron 
sea, then the chemical potential becomes locally shifted so that p ^ p + ef(t). So long 
as the field 0(t) is varied at a rate slowly compared with the Kondo temperature, the 
constraint field will always track with the chemical potential, and since the constraint 
field is pinned to the chemical potential, A ^ A + efit). In the process, the constraint term 
will become 


m/U) -Q)^ miU) -Q) + efmfU) - Q). 


(50) 
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Sinee the f-eleetrons now eouple to the external potential e4> we have to aseribe a physieal 
charge e = -\e\ to them. By contrast, the -Q term in the constraint must be interpreted 
as a “background positive charge” \e\Q = \e\ per site. These lines of reasoning indicate 
that we should think of the Kondo effect as an many-body ionization phenomenon in 
which the neutral local moment splits up into a negatively charged heavy electron and a 
stationary positive background charge we can associate with the formation of a Kondo 
singlet. 


3.4 The Composite nature of the f-electron 

The matrix Green’s function of the Kondo lattice reminds us of the Nambu Green’s function in 
superconductivity. It is given by 


^k(T) = -<iAko-(T)<A'^ko-(0)) 


Gc(k, r) Gc/(k, r) 
G/c(k,T) G/(k,T) 


(51) 


where Gc(k,T) = -(ck(T)c'*’ko-(0)), Gc/(k,r) = -(ck(T)/'*'ko-(T)) and so on. The anomalous off- 
diagonal members of this Green’s function remind us of the Gor’kov functions in BCS theory, 


and develop with the coherent hybridization. Using the two component notation (35), this 
Green’s function can be written 


F.T. 

Q^{r) = -{dr + h^r' - 

where F.T denotes a Fourier transform in imaginary time {dr 

&k{z) = {z- 


^ 0k{^^n} — ) (52) 

-i(On), or more explicitly. 


1 _ 

z - 6k -y ' 

-1 

'G,(k,z) 

G,/(k,z)' 


,-y z-a) 


^G/c(k,z) 

Gf(k,z), 


1 


(z-ek)(z-d)-y2 


z- A 
V 


V 

z-ek 


(53) 


where we have taken the liberty of analytically extending icor z into the complex plane. 
Now we can read off the Green’s functions. In particular, the “hybridized” conduction electron 
Green’s function is 

z- A 


G,(k,z) 


1 


(z-ek)(z-A)-V^ 

1 


z-e,-i:Az) 


(54) 


which we can interpret physically as conduction electrons scattering off resonant f-states at each 
site, giving rise to a momentum-conserving self energy 


i:c(z) 


V 


V 


y 2 


0 ( 1 ) 


z- A 


(55) 


We see that the Kondo effect has injected a resonant scattering pole at energy z = T in the 
conduction electron self-energy. This resonant scattering lies at the heart of the Kondo effect. 
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3.4.1 An absurd digression: the nuclear Kondo effect 

The appearanee of this pole in the seattering raises a vexing question in the Kondo effeet: what is 
the meaning of the f-electron ? This might seem like a dumb question, for in eleetronic materials 
the Kondo effect certainly involves localized f electrons, and surely, we can interpret this pole 
as as the adiabatic renormalization of a hybridized band-structure. This is certainly true. Yet 
as purists, we do have to confess that our starting model, was a pure Kondo lattice model with 
only spin degrees of freedom: they could even have been nuclear spins! 

This might seem absurd, yet nuclear spins do couple antiferromagnetic ally with conduction 
electrons to produce nuclear antiferromagnetism. Leaving aside practical issues of magnitude, 
we can learn something from the thought experiment in which the the nuclear spin coupling to 
electrons is strong enough to overcome the nuclear magnetism. In this case, resonant bound- 
states would form with the nuclear spin lattice giving rise to charged heavy electrons, presum¬ 
ably with an expanded Fermi surface. 

From this line of argument we see that while it’s tempting to associate the heavy fermion with a 
physical f- or d- electron localized inside the local moment, from a renormalization group per¬ 
spective, the heavy electron is an emergent excitation: a fermionic bound-state formed between 
the conduction sea and the neutral localized moments. This alternate point-of-view is useful, 
because it allows us to contemplate the possibility of new kinds of Kondo effect into states that 
are not adiabatically accessible from a band insulator or metal. 


3.5 Cooper pair analogy 


There is a nice analogy with superconductivity which helps to understand the composite nature 
of the heavy electron. In a superconductor, electron pairs behave as loose composite bosons 
described by the relation 

= -Fix -x'). (56) 


Here Fix - x') = -iT ^^(l)i/fj^(2)) is the anomalous Gor’kov Greens function which determines 
the Cooper pair wavefunction, extended over the coherence length ^ ~ Vp/Tc. A similar phe¬ 
nomenon takes place in the Kondo effect, but here the bound-state develops between spins and 
electrons, forming a fermion, rather than a boson. For the Kondo lattice, its perhaps more use¬ 
ful to think in terms of a screening time tk ~ hlTp, rather than a length. Both the Cooper pair 
and heavy electron involve electrons that span decades of energy up to a cutoff, belt the Debye 
energy cu/) in superconductivity or the (much larger) bandwidth D in the Kondo effect [ 34}^ . 
To follow this analogy in greater depth, recall that in the path integral the Kondo interaction 
factorizes as 

v(cV.) + (/V«)v + A^, (57) 


L t c 


J 


SO by comparing the right and left hand side, we see that the composite operators SpaCp and 
fjS ap behave as a single fermion denoted by the contractions: 
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1 ^ 

(y)A, (58) 

' ' 




Composite Fermion 


Physieally, this means that the spins bind high energy eleetrons, transforming themselves into 
eomposites whieh then hybridize with the conduetion eleetrons. The resulting “heavy fermions” 
can be thought of as moments ionized in the magnetically polar electron fluid to form mobile, 
negatively charged heavy electrons while leaving behind a positively charged “Kondo singlet”. 
Microscopically, the many body amplitude to scatter an electron off a local moment develops a 
bound-state pole, which for large N we can denote by the diagrams: 



0(1 m 


The leading diagram describes a kind of “condensation” of the hybridization field; the second 
and higher terms describe the smaller 0(1 /N) fluctuations around the mean-field theory. 

By analogy with superconductivity, we can associate a wavefunction associated with the tem¬ 
poral correlations between spin-flips and conduction electrons, as follows 

^ ^ C/3(r)^flair') = g(T - T')faiT'). (59) 

where the spin-flip correlation function g{T-T') is an analogue of the Gor’kov function, extend¬ 
ing over a coherence time Tk ~ ^ITk- Notice that in contrast to the Cooper pair, this composite 
object is a fermion and thus requires a distinct operator fa for its expression. 



4 Heavy Fermion Superconductivity 


We now take a brief look at heavy fermion superconductivity. There are a wide variety of 
heavy electron superconductors, almost all of which are nodal superconductors, in which the 
pairing force derives from the interplay of magnetism and electron motion. In the heavy 
fermion compounds, as in many other strongly correlated electron systems superconductiv¬ 
ity frequently develops at the border of magnetism, near the quantum critical point where the 
magnetic transition temperature has been suppressed to zero. In some of them, such as UPta 


(Tc=0.5K) p6| the superconductivity develops out of a well-developed heavy Fermi liquid, 
and in these cases, we can consider the superconductor to be paired by magnetic fluctuations 
within a well-formed heavy Fermi liquid. However, in many other superconductors, such as 
UBei 3 (Tc=lK) [[^[^, the 115 superconductors CeCoInj (Tc=2.3K) [ [39| , CeRhlnj under 
pressure (Tc=2K) fT^, NpAl2Pd5(T,=4.5K) ig and PuCoGag (Tc=18.5K) the su¬ 

perconducting transition temperature is comparable with the Kondo temperature. In many of 
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these materials, the entropy of eondensation 



(60) 


ean be as large as (1/3)7? In 2 per rare earth ion, indicating that the spin is, in some-way entan¬ 
gling with the conduction electrons to build the condensate. In this situation, we need to be able 
to consider the Kondo effect and superconductivity on an equal footing. 



CeColns 



Fig. 12: (a) Phase diagram of 115 compounds CeMln^, adapted from /[43]/, showing magnetic 
and superconducting phases as a function of alloy concentration, (b) Sketch of specific heat co¬ 
efficient of CeColn^, (with nuclear Schottky contribution subtracted), showing the large entropy 
of condensation associated with the superconducting state. (After Petrovic et al 2001 /|39]/). 


4.1 Symplectic spins and SP (N). 

Although the SU(N) large N expansion provides a very useful description of the normal state 
of heavy fermion metals and Kondo insulators, there is strangely, no superconducting solution. 
This short-coming lies in the very structure of the S U(N) group. S U(N) is perfectly tailored 
to particle physics, where the physical excitations - the mesons and baryons appear as color 
singlets, with the meson a a qq quark-antiquark singlet while the baryon is an A-quark singlet 
q\q 2 - ■ -qN, (where of course A = 3 in reality). In electronic condensed matter, the meson 
becomes a particle-hole pair, but there are no two-particle singlets inSU (A) beyond A = 2. The 
origin of this failure can be traced back to the absence of a consistent definition of time-reversal 
symmetry in SU (A) for A > 2. This means that singlet Cooper pairs and superconductivity can 
not develop at the large A limit. 


A solution to this problem which grew out an approach developed by Read and Sachdev [441 for 


frustrated magnetism, is to use the symplectic group S P(N), where A must be an even number 
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1^1^. This little-known group is a subgroup of S U(N). In faet for N = 2, SU(2) = S P{2) 
are identical, but they diverge for higher N. For example, S U(4) has 15 generators, but its 
symplectic sub-group S P{4) has only 10. At large N, S P{N) has approximately half the number 
of generators of SU(N). The symplectic property of the group allows it to consistently treat 
time-reversal symmetry of spins and it also allows the formation of two-particle singlets for any 
N. 


One of the interesting aspects of S P(N) spin operators, is their relationship to pair operators. 
Consider SP(2) = S U{2): the pair operator is since this operator is a sin¬ 
glet, it commutes with the spin operators, S] = = 0 which, since W and are 

the generators of particle-hole transformations, implies that the S U (2) spin operator is particle- 
hole symmetric. It is this feature that is preserved by the S P{N) group, all the way out to 
A ^ oo. In fact, we can use this fact to write down an SP{N) spins as follows: an S U{N) 
spin is given by = f'^afp- Under a particle hole transformation /„ ^ Sgn(cr)/_Q,. If we 

take the particle-hold transform of the S U(N) spin and add it to itself we obtain an S P(N) spin. 


Sap = fafp + Sgn(Q:yS)/_^/^_„, 

(61) 


Symplectic Spin operator 


where the values of the spin indices are a,/? 6 {+1/2,..., +N/2}. This spin operator commutes 
with the three isospin variables 


T 3 = rif - A/2, ^ 

a>0 a>0 


(62) 


With these local symmetries, the spin is continuous invariant under SU (2) particle-hole rota¬ 
tions fa ufa + vSgnQ'/'*'_Q,, whcrc \u^\ + |v^| = 1, as you can verify. To define an irreducible 
representation of the spin, we also have to impose a constraint on the Hilbert space, which in its 
simplest form is T 3 = r* = 0, equivalent to 2 = A/2 in the S U(N) approach. In other-words, 
the s-wave part of the f-pairing must vanish identically. 



T 


7. 



Andrei and Coleman 


1989 


Fig. 13: Phase diagram for the two-dimensional Kondo Heisenberg model, derived in the 
S P(N) large A approach, adapted from /|?7|/, courtesy Rebecca Flint. 
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4.2 Superconductivity in the Kondo Heisenberg Model 


Let us take a look at the way this works in a nearest neighbor “Kondo Heisenberg model” [471, 


H = H, + Hk + Hm. 


(63) 


Here He = 2ko- e^kC^o-Cko- deseribes the eonduetion sea, whereas Hk and Hm are the Kondo and 
Heisenberg (RKKY) interaetions, respectively. These take the form 

Hk = ^ ^ C^jaCjpSpc^U) ^ ^ jaf j-a){fj-l3Cjp)) 

j iJ 

Hm = ^ Xi (64) 

('.;■) j 


where we’ve introduced the notation a = Sgn(a;) and have shown how the interactions are 
expanded in into particle-hole and particle-particle channels. Notice how the interactions are 
equally divided between particle-hole and particle-particle channels. When we carry out the 
Hubbard Stratonovich decoupling, in each of these terms, we obtain 


Hk 

Hh 


ny .|2 + 

E + + H.c]+W ' 


y) \tijfiofja + Aij&fiafi-o + H.cj + N 

(ij) 


Jk , 
\tij? + M,-/ 


J, 


H 


(65) 


At each site, we can always rotate the f-electrons in particle-hole space to remove the “Kondo 
pairing” component and set = 0, but the pairing terms in the Heisenberg component can 
not be eliminated. This mean-field theory describes a kind of Kondo stabilized spin-liquid [ |47| . 
The physical picture is as follows: in practice, a spin-liquid is unstable to magnetism, but its 
happy co-existence with the Kondo effect brings its energy below that of the antiferromagnet. 
The hybridization of the / with the conduction sea converts the spinons of the spin-liquid into 
charged fermions. The terms describe various kind of exotic density waves. The zl,y terms 
now describe pairing amongst the composite fermions. 

To develop a simple theory of the superconducting state, we restrict our attention to uniform, 
static saddle points, dropping the tij. Lets look at the resulting mean-field theory. In two dimen¬ 
sions, this becomes 


H= 


k,Q;>0 


fk^a Vti 

Vti w-t + Ah^^ti 


^ki 


Aq 


J \jKaJ 


|V|2 . JA 


+ ^ + 2 


Hi 
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J, 


H 


( 66 ) 


where 


^ ka ~ ka^ (^^-k,-a\ f ka kaj f—k—a) 


(67) 


are Nambu spinors for the conduction and f-electrons. The vector l^of Lagrange multipliers 
couples to the isospin of the f-electrons: stationarity of the Free energy with respect to this 
variable imposes the mean-field constraint that iprf) = 0. The function Ank = Ak(coskji - 
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eos ky) is the f-electron pair wavefunetion. Here we’ve ehosen a d-wave form-faetor. For this 
choiee, the loeal f pair density automatically vanishes and so we need only choose w = (0,0, A), 
where A couples to T3 (imposing the constraint ny = N/2). We could have also tried an extended 
extended s-wave pair wavefunetion, but in this case, the induced s-wave pair density becomes 
finite, and the effect of the w constraint is to suppress the transition temperature. By seeking 
stationary points in the free energy with respect to variations in Ah, V and A one can derive 
the phase diagram for d-wave pairing, shown in Fig. 13 The mean-field theory shows that 
superconductivity develops at the interface between the Fermi liquid and the spin liquid. 


5 Topological Kondo Insulators 


One of the areas of fascinating development in the last few years, is the discovery that Kondo 
insulators can develop topological order to form a Topological Kondo insulator Topological 
order refers to the idea that a quantum mechanical ground-state can develop a non-trivial topol¬ 
ogy. One of the defining features of topological ground-states is the development of protected 
surface states. The best known example of topological order the integer quantum Hall effect, 
where an integer filled Landau level develops topological order that is is responsible for the 
robust quantization of the Quantum Hall effect [[4^|-(50|. In a remarkable series of discoveries 


in 2006, [51-581 it became clear that strong spin orbit coupling can play the role of a synthetic 
magnetic field, so that band insulators can also develop a non-trivial topology while preserving 
time-reversal symmetry. Such Z2 topological band insulators are defined by a single topolog¬ 
ical Z2 = +1 index that is positive in conventional insulators, but reverses in topological Z2 
insulators. This topological feature manifests itself through the formation of robust conducting 
surface states. 

In 2007, Liang Fu and Charles Kane showed that if an insulator has both time reversal and 


inversion symmetry [^, this Z2 index is uniquely determined by the the parities of the 
Bloch states at the high symmetry points T,- of the valence band which determine a “Z2 index” 


Z2 




-I -1 conventional insulator 
-1 topological insulator 


( 68 ) 


Fu Kane formula for the Z2 index of topological insulators 


where d(r,) = dm is the the product of the parities of the of the occupied bands at the 
high-symmetry points in the Brillouin zone. This formula allows one to determine whether an 
insulator state is topological, merely by checking whether the index Z2 = -1, without a detailed 
knowledge of the ground-state wavefunetion. 

It used to be thought that Kondo insulators could be regarded as “renormalized silicon”. The 
discovery of topological insulators forced a re-evaluation of this viewpoint. The large spin orbit 
coupling, and the odd-parity of the f-states led to the proposal, by Dzero, Sun, Galitski and the 


author, [59 j that Kondo insulators can become topologically ordered. The Fu-Kane formula has 
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a speeial signifieanee for Kondo insulators, whieh eontain odd parity f-eleetrons hybridizing 
with even parity d-electrons. Each time an f-electron crosses through the band-gap, exchanging 
with a conduction d-state, this changes the Z 2 index, making it highly likely that certain Kondo 
insulators are topological. The oldest known Kondo insulator SmBg, discovered almost 50 years 


ago was well known to possess a mysterious low temperature conductivity plateau [ |6Q|[6T| , and 
the idea that this system might be a topological Kondo insulator provided an exciting way 
of explaining this old mystery. The recent observation of robust [^63| conducting surface 
states in the oldest Kondo insulator SmBgSupports one of the key elements of this prediction, 
prompting a revival of interest in Kondo insulators as a new route for studying the interplay of 
strong interactions and topological order. 

SmBg is really a mixed valent system, which takes us a little beyond the scope of this lec¬ 
ture. One of the other issues with SmBg, is that its local crystal field configuration is likely 


to be a Eg quartet state [64|, rather than a Kramers doublet. Nevertheless, key elements of its 
putative topological Kondo insulating state are nicely illustrated by a spin-orbit eoupled Kondo- 
Heisenberg model, describing the interaction of Kramer’s doublet f-states with a d-band. The 
model is essentially identical with (Eq. 


(69) 


kcr^kcr dx ^ ^ ja^ Jh ^ ^ 

kcr ;■ 1, / 


with an important modification that takes into account the large spin-orbit coupling and the 
odd-parity of the f-states. This forces the local Wannier states Wja that exchange spin with the 
local moment to be odd parity combinations of nearest neigbour conduction electrons, given by 

^ io-<PtTa(Ri - Ry) (70) 

i,(T 


We’ll consider a simplified model with the form factor 


0(R) = 


-iR ‘I, Re n.n 
0 otherwise 


(71) 


This form factor describes the spin-orbit mixing between states with orbital angular momentum 
I differing by one, such as / and d or p and s orbitals. The odd-parity of the form-factor 
0(R) = -0(-R) derives from the odd-parity /- orbitals, while the prefactor -i ensures that 
the hybridization is invariant under time-reversal. The Fourier transform of this Form factor, 
^(k) = 2 r ^(R)e''^ is then 

0(k) = 4 • c? (72) 

where the s-vector 4 = (sinki,sink 2 >sink 3 ) is the periodic equivalent of the unit momentum 
vector k. Notice how 5(r,) = 0 vanishes at the high symmetry points. 

The resulting mean-field Hamiltonian takes the form 


Htki - E i//\h(k)if/k + Ns 

k 


Y + -r-^Q 

Jk Jh 


(73) 
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where iIj\ = {c\a,f\a-) and 


/z(k) = 


/ 


fk 


ya ■ 4 




(74) 


while f/k = 2tf{Cx+Cy + c^ + X (c/ = eos kt) is the dispersion of the f-state resulting from a mean- 
field deeoupling of the intersite Heisenberg coupling in the particle-hole channel. For small k, 
the hybridization in Hamiltonian fi(k) takes the form Vc? • k, a form which closely resembles 
the topologically non-trivial triplet p-wave gap structure of superfluid He-3B. Like He-3B, the 
hybridization only develops at low temperatures, making SmBg an adaptive insulator. 


a) 

Z2(a)=+1 



b) 


Z2(b)=(-ip 



Fig. 14: (a) When the d-band is above the filled f-band, a trivial insulator is formed, (b) When 
the d-band crosses the f-band at the three X-points, the Z 2 parity changes sign, giving rise to a 
topological insulator. 


Let us for the moment treat fi(k) as a rigid band structure. Suppose the f-band were initially 
completely filled, with a completely empty d-band above it. (See Fig. [T^ a). This situation 
corresponds to a conventional band insulator with Zj = Next, let us lower the d-conduction 
band until the two bands cross at a high symmetry point, causing the gap to close, and then 
to re-open. We know, from dHvA studies of the iso-electronic material LaBg [651 (whose 
band-structure is identical to SmBe but lacks the magnetic f-electrons), from ARPES studies 
[?, 66-681, that in SmBg, the d-band crosses through the Fermi surface at at the three X points. 
Once the d-band is lowered through the f-band around the three X points, the odd-parity f-states 
at the X point move up into the conduction band, to be replaced by even-parity d-states. This 
changes the sign of Z 2 ^ (-1)^ = -1, producing a topological ground-state. Moreover, since 
there are three crossing, we expect there to be three spin-polarized surface Dirac cones. 

We end by noting that at the time of writing, our understanding of the physics SmBg is in 
rapid flux on both the experimental and theoretical front. Spin resolved ARPES [ [69| mea¬ 
surements have detected the presence of spin-textures in the surface Fermi surfaces around the 
surface X point, a strong sign of topologically protected surface states. Two recent theoreti- 
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eal works [70,711 have shown that the spin textures seen in these experiments are eonsistent 
with a spin-quartet ground-state in SmBe- Despite this progress, consensus on the topological 
nature of SmBg has not yet been achieved, and competing groups have offered alternate inter¬ 


pretations of the data, including the possibility of polarity-driven surface metallicity [721 and 


and Rashba-split surface states, both, of a non-topological origin. [73[. Another area of exper¬ 
imental controversy concerns the possible de-Haas van Alphen oscillations created by surface 
topological excitations, with one report of the detection of surface de Haas van Alphen sig¬ 


nals [74[ and a recent, very remarkable report of bulk de Haas van Alphen signals associated 


with unhybridized, quantum critical d-electrons [75 j. 


6 Co-existing magnetism and Kondo effect 

In this short lecture. I’ve given a quick introduction to the paramagnetic phases of heavy fermion 
systems. One of the of major open questions in heavy fermion and Kondo lattice physics con¬ 
cerns the physics of magnetism, and the right way to describe the development of magnetism 
within these materials. There is growing evidence that magnetism and the Kondo effect can co¬ 
exist, sometimes homogeneously, and sometimes inhomogeneously. For example. In the 115 
superconductor CeRhlnj there is evidence for a microscopic and homogeneous coexistence of 
local moment magnetism and heavy fermion superconductivity under pressure [ |7^ ; By con¬ 
trast, in the geometrically frustrated CePdAl [77}78|, two thirds of the Cerium sites sponta¬ 


neously develop magnetism, leaving the other third to undergo a Kondo effect [79[. What is the 
right way to describe these co-existent states? 

One possibility that 1 have worked on with Aline Ramires [ [80|[M] [ is the use of a supersymmetric 
is the use of a “supersymmetric” spin representation of the spin 


S a/i - f\fp + b\bp 


(75) 


where the f\ and „ are fermionic and bosonic creation operators. Such a representation per¬ 
mits in principle, the existence of “two fluid” ground-states, involving a Gutzwiller projection 
of bosonic and fermionic wavefunctions 


\W) = Pg\Wp)\Wb), 


(76) 


where \ Wp) is the fermionic component of the wavefunction describing the Kondo quenched 
local moments while | Wp) describes the formation of long-range magnetic correlations within a 
bosonic RVB wavefunction, while 


=/n 


^ iejiriB+rif-l) 

2n 


(77) 


is a Gutzwiller projection operator onto the state with one spin per site. We have been trying 
to describe such mixed state wavefunctions in the large-A limit, seeking saddle point solutions 


where a bosonic and fermionic fluid co-exist [81 j. One of the ideas that emerges from this kind 
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of approaeh, is the possibility that the soft modes at a Quantum Critieal point might develop 


fermionie character, a kind of emergent supersymmetry [82|. 
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